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Signalauswerteverfahren zur Detektion von QRS-Komplexen in Elek- 
trokardiogramm-Signalen 

Die Erfindung betrifift ein Signalauswerteverfahren zur Detektion von QRS- 
5 Komplexen in Elektrokardiogramm-(EKG- und IEGM-)Signalen. 

Zum Hintergrund der Erfindung ist festzuhalten, daB die automatische 
Analyse von EKG-Signalen zur Perfektionierung der Funktionsfahigkeit 
^ von Herzschrittmachern und Defibrillatoren eine immer groBere Rolle 
10 spielt. Solche kardiologischen Gerateimplantate jungerer Bauart bieten 
dementsprechend auch die Fahigkeit einer EKG-Analyse. In diesem Zu- 
sammenhang hat die Erkennung von QRS-Komplexen und R-Zacken in 
EKG-Signalen einen uberragenden Stellenwert. Diese Bedeutung resultiert 
aus der vielfaltigen Verwendbarkeit der Information uber das zeitliche 
15 Auftreten des QRS-Komplexes, z. B. bei der Untersuchung der Herzraten- 
variabilitat, bei der Klassifikation und bei der Datenkompression sowie als 
Basissignal fur sekundare Anwendungen. Oberhaupt nicht oder falsch de- 
tektierte QRS-Komplexe und R-Zacken sind dabei problematisch im Hin- 
@ ^ blick auf die Leistungsfahigkeit der auf die Detektion folgenden Verarbei- 
20 tungs- und Analysestufen. 

Einen breiten Oberblick uber bekannte Signalauswerteverfahren zur Detek- 
tion von QRS-Komplexen in EKG-Signalen bietet der Fachaufsatz von 
Friesen et al. "A Comparison of the Noise Sensitivity on Nine QRS De- 
25 tection Algorithms" in IEEE Transaction on Biomedical Engineering, Vol. 
37, Nr. 1, Januar 1990, Seiten 85 - 98. Die dort vorgestellten Signalaus- 
werte-Algorithmen basieren durchgehend auf einer Auswertung der 
Amplitude, der ersten Ableitung des Signals sowie dessen zweiter Ablei- 
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tung. Bei den vorgestellten Algorithmen wird unterschieden zwischen sol- 
chen, die eine Analyse der Amplitude und der ersten Ableitung vornehmen, 
solchen, die nur die erste Ableitung bewerten und solchen, die die erste und 
zweite Ableitung in Betracht ziehen. Kurz umrissen wird bei alien Algo- 
5 rithmen abgepruft, ob der jeweilige Signalparameter vorgegebene Schwel- 
len uber- oder unterschreitet, wonach bei Auftreten eines solchen Ereignis- 
ses mit einem vorgegebenen Schema das Auftreten weiterer definierter Er- 
eignisse abgepruft und bei Erfullung bestimmter Kriterien auf das Vorlie- 
, gen eines QRS-Komplexes geschlossen wird. 

10 

Ein anderer Aspekt bei der Signalauswertung zur Detektion von QRS- 
Komplexen ist bei der Implementierung solcher Verfahren in implantierten 
Herzschrittmachern zu beachten. Im Hinblick auf die naturgemafi bei sol- 
chen Geraten gegebenen Beschrankungen bei der Energieversorgung und 
15 der Rechenkapazitat ist darauf zu achten, daB zur Detektion von QRS- 
Komplexen moglichst einfache Algorithmen mit moglichst wenigen Re- 
chenoperationen auf der Basis von ganzen Zahlen anstatt von reellen Zah- 
len durchzufuhren sind. 

" 20 Ferner wurden zur QRS-Erkennung Signalverarbeitungsmeuioden aus den 
Gebieten der linearen und nicht-linearen Filterung, der Wavelet- 
Transformation, der kunstlichen neuronalen Netze und der genetischen Al- 
gorithmen angewendet. Bei hohen Signal-Stor-Abstanden und nicht patho- 
logischen Signalen, also im Falle guter Signalverhaltnisse, erzielten diese 
25 Auswerteverfahren zuverlassige Ergebnisse. Waren keine solchen Verhalt- 
nisse gegeben, konnte die Leistungsfahigkeit der Auswerteverfahren unter 
Umstanden drastisch absinken, was im Hinblick auf eine zuverlassige 
Funktionsweise eines Schrittmachers natiirlich nicht akzeptabel ist. 



SchlieBlich ist aus der alteren, jedoch nachveroffentlichten deutschen Pa- 
tentanmeldung 100 11 733.3 der Anmelderin eine QRS-Erkennung auf der 
Basis von Nulldurchgangszahlungen bekannt. Dabei werden folgende Ver- 
fahrensschritte durchgefiihrt: 

- Abtasten des Signals und Umwandlung in diskrete, zeitlich aufein- 
anderfolgende Signalwerte, 

- HochpaBfilterung der abgetasteten Signalwerte, 

- Bestimmung des Vorzeichens jedes Signalwerts, 

- laufendes Abpriifen der Vorzeichen aufeinanderfolgender Signal- 
werte auf Vorliegen eines Nulldurchgangs zwischen zwei aufeinan- 
derfolgenden Signalwerten, 

- Ermittlung der Anzahl der Nulldurchgange in einem definierten 
Segment der aufeinanderfolgenden Signalwerte, und 

- Vergleichen der ermittelten Nulldurchgangszahl mit einem definier- 
ten Schwellenwert, wobei ein Unterschreiten des Schwellenwertes 
signifikant fur das Vorliegen eines QRS-Komplexes in dem defi- 
nierten Segment des Signalverlaufes ist. 

Um dabei im Bereich aufierhalb der QRS-Komplexe eine signifikant hohe 
Nulldurchgangszahl zu erhalten, wird dort zu dem Hoch- oder BandpaB- 
gefilterten und quadrierten EKG-Signal ein hochfrequentes Uberlagerungs- 
signal mit im Vergleich zur Amplitude des QRS-Komplexes niedriger 
Amplitude addiert. 

Diese Vorgehensweise steht der eingangs erlauterten Forderung nach mog- 
lichst einfachen Algorithmen natiirlich entgegen. 



Ausgehend von den geschilderten Problematiken liegt der Erfindung die 
Aufgabe zugrunde, ein Signalauswerteverfahren zur Detektion von QRS- 
Komplexen in EKG-Signalen anzugeben, das mit vergleichsweise geringer 
Rechenkapazitat und auch bei problematischen Signalverhaltnissen einsetz- 
bar ist, dabei jedoch zuverlassige Detektionsergebnisse liefert. 

Diese Aufgabe wird durch die im Patentanspruch 1 angegebenen kenn- 

zeichnenden Verfahrensschritte wie folgt gelost: 

- Abtasten des EKG-Signals und Umwandlung in diskrete, zeitlich 

aufeinanderfolgende Signalwerte, 

- Vergleichen der Signalwerte mit einer adaptiv aus diesen ermittelten 
SchweUenwertfunktion, 

- Feststellen einer Haufigkeitszahl innerhalb eines definierten Seg- 
mentes der aufeinanderfolgenden Signalwerte, mit der vorzugsweise 
die Betrage der Signalwerte die SchweUenwertfunktion unterschrei- 
ten, und 

- Vergleichen der ermittelten Haufigkeitszahl mit einem definierten 
Grenzwert, wobei ein Unterschreiten des Grenzwerts signifikant fur 
das Vorliegen eines QRS-Komplexes in dem definierten Segment 
des EKG-Signals ist. 

Kernpunkt des erfmdungsgemaBen Verfahrens ist die Anwendung eines 
Schwellenwertvergleichs und einer nachfolgenden Haufigkeitszahlung, 
wobei die Morphologie des QRS-Komplexes ausgenutzt wird. Der QRS- 
Komplex im EKG-Signal ist namlich durch eine relativ hochamplitudige 
Oszillation charakterisiert, die den Signalverlauf ausgepragt von der in der 
Regel verrauschten und mit einem Offset beaufschlagten Basis-Linie des 
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Elektrokardiogramms wegfuhrt. Die Frequenz dieser kurzen Oszillation 
liegt in einem Bereich, in dem andere Signalkomponenten, wie die P- und 
T- Welle, nur geringen EinfluB ausiiben und vorzugsweise durch Vorfilte- 
rung - z. B. eine Hoch- oder BandpaBfilterung - entfernt werden konnen. 

5 Nach der Unterdriickung dieser tieffrequenten Signalkomponenten ergeben 
sich aufgnind von hoherfrequentem Rauschen Signalschwankungen urn die 
Basis-Linie, die in dem Bereich, wo kein QRS-Komplex auftritt, das EKG- 
Signal dominieren. Der QRS-Komplex erscheint dann in diesem Signal- 
kontext als langsame, nur kurz andauernde Schwingung hoher Amplitude, 
1 10 die von der Basis-Linie des EKG-Signals signifikant wegfuhrt. Werden also 
die Signalwerte mit einer das Signalrauschen reprasentierenden Schwel- 
lenwertfunktion verglichen, so sind die Betrage der Signalwerte auBerhalb 
des QRS-Komplexes meist kleiner als diese Schwellenwertfunktion. Inso- 
weit ist die Haufigkeitszahl, mit der die Betrage der Signalwerte die 

15 Schwellenwertfunktion unterschreiten, grofi. Im Bereich des QRS- 
Komplexes uberschreiten mm die Betrage der Signalwerte die Schwellen- 
wertfunktion, und zwar signifikant. Es stellt sich also eine geringe Haufig- 
keitszahl fur Unterschreitungen der Schwelle im Zuge des QRS-Komplexes 
ein. Durch ein Vergleichen der ermittelten Haufigkeitszahl mit einem defi- 

20 nierten Grenzwert kann also der QRS-Komplex selektiert werden. Ein Un- 
terschreiten des Grenzwertes signalisiert namlich das ftir den QRS- 
Komplex typische Oberschreiten der Schwellenwertfunktion. 



Das erfindungsgemaBe Verfahren der QRS-Komplex-Erkennung hat sich 
25 als robust gegenuber Rauschstorungen imd rechentechnisch besonders ein- 
fach realisierbar erwiesen. Insoweit eignet es sich besonders fur die Imple- 
mentierung bei der Realzeit- Analyse von EKG-Signal-Morphologien in 
Herzschrittmachern. 
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Die bereits erwahnte HochpaBfilterung wird vorzugsweise mit einer unte- 
ren Durchlafifrequenz von 18 Hz vorgenommen. Damit k6nnen die tieffre- 
quenten Komponenten, wie die P- und T-Welle sowie eine Grundliniendrift 
5 unterdruckt werden. Ferner wird der QRS-Komplex dadurch zur Signal- 
komponente mit der niedrigsten Frequenz, die das Signal wShrend ihres 
Auftretens dominiert. 

Zur VergroBerung des Signal-Stor-Abstandes kann ferner vorgesehen sein, 
10 die Signalwerte vor dem Vergleichen mit der Schwellenwertfunktion und 
der Haufigkeitszahl zu quadrieren. Damit werden kleine Signalbetrage re- 
lativ zu grofieren Signalbetragen abgeschwacht, was die Detektierbarkeit 
des QRS-Komplexes weiter verbessert. 

15 Der Wert der SchweUenwertfunktion wird vorzugsweise adaptiv aus einer 
flieBenden Mittelung der BandpaB-gefilterten und quadrierten Signalbetra- 
ge bestimmt. 

Im folgenden wird das erfindungsgemaBe Verfabren anhand der beigefug- 
20 ten Zeichnungen in einem Ausfuhrungsbeispiel naher erlautert. Es zeigen: 

Fig. 1 eine hochst schematische Darstellung des Signalverlaufs eines 
QRS-Komplexes in einem EKG-Signal, und 

25 Fig. 2 ein Strukturschema des erfmdungsgemaBen Signalauswertever- 
fahrens zur Detektion von QRS-Komplexen in EKG-Signalen. 



Wie aus Fig. 1 deutlich wird, besteht ein schematisch dargestellter QRS- 



Komplex aus einer relativ hochamplitudigen Oszillation urn die NuU-Linie 
1, die das EKG-Signal 4 in der Q-Zacke 6 zuerst in negativer Richtung iiber 
die Null-Linie 1 wegfuhrt. AnschlieBend wird das EKG-Signal 4 in der R- 
Zacke 5 mit einem steilen Anstieg in den positiven Bereich und einem an- 
5 schlieCenden steilen Abfall wiederum in den negativen Bereich unter Aus- 
bildung der S-Zacke 7 ubergefuhrt. 

Real ist das EKG-Signal 4 mit einem bestimmten Pegel verrauscht, was in 
Fig. 1 durch den strichlierten Signalverlauf angedeutet ist. Wird nun dieses 
^ 10 verrauschte Signal 4 abgetastet und in diskrete, zeitlich aufeinanderfolgen- 
de Signalwerte umgewandelt und BandpaB-gefiltert, so konnen diese 
Signalwerte mit der in Fig. 1 schematisch als mit Kreuzchen unterlegte Li- 
nie dargestellten SchweUenwertfunktion K(n) verglichen werden. Wie sich 
mm aus Fig. 1 anschaulich und modellhaft ableiten laBt, liegt im Bereich 
15 aufierhalb des QRS-Komplexes der Betrag des EKG-Signals meist unter- 
halb dieser SchweUenwertfunktion K(n). Es ergeben sich also beispielswei- 
se im Bereich Nl signifikant hohe Haufigkeiten fur Signalbetrage I x(n) | 
imterhalb der SchweUenwertfunktion K(n). 

^20 Im Bereich der Q-Zacke 6 weicht nun der Betrag des EKG-Signals sehr 
stark von der SchweUenwertfunktion K(n) in positive Richtung ab. Die 
Haufigkeitszahl D(n) innerhalb des Segmentes N2 des QRS-Komplexes 6 
fur dieses Ereignis ist also wesentlich geringer als die Haufigkeitszahl in- 
nerhalb des Segments Nl. Insoweit kann die Haufigkeitszahl D(n) fur die 
25 Erkennung des QRS-Komplexes herangezogen werden. Dieser wird dann 
als vorliegend erkannt, wenn die Haufigkeitsanzahl D(n) einen defmierten 
Grenzwert 0 unterschreitet. 



• 
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Als ein Kernpunkt der Erfmdung gegentiber dem Stand der Technik ist an 
dieser Stelle hervorzuheben, daB aufgrund der Feststellung der erSrterten 
HSufigkeitszahl und deren Vergleich mit einem definierten Grenzwert nicht 
- wie im Stand der Technik - die Amplitude des EKG-Signals daraufhin 

5 uberpruft wird, ob ein bestimmter Schwellenwert absolut uberschritten 
wird, urn daraus auf den QRS-Komplex zu schlieBen. Vielmehr wird ge- 
wissermaBen iiberpruft, wie lange das EKG-Signal eindeiitig auf einer fur 
das Vorliegen eines QRS-Komplexes sprechenden Seite einer Schwellen- 
wertfunktion „verharrt". Erst das Vorliegen einer bestimmten Zeitdauer 

10 dieses Zustandes wird zum RiickschluB auf das Vorliegen eines QRS- 
Komplexes herangezogen. Demzufolge werden nur kurze, starke MeB- 
schwankungen nicht als (falsche) QRS-Komplexe (sogenannte „False posi- 
tive"-Fehler) erkannt. 

1 5 Der detaillierte Ablauf des erfindungsgemaBen Auswerteverfahrens wird 
nun im folgenden anhand der Fig. 2 naher erlautert. 

Das EKG-Signal 4 wird abgetastet und in diskrete, zeithch aufeinanderfol- 
gende Signalwerte x(n) umgewandelt. Die Abtastrate betragt beispielsweise 

20 f T = 360 Hz, d. h., daB das EKG-Signal in eine Folge von 360 MeBwerten 
pro Sekunde umgewandelt wird. Das abgetastete EKG-Signal x(n) wird 
dann eingangsseitig einer BandpaBfilterung BP unterzogen, die der Entfer- 
nung aller Signalkomponenten, die nicht zum QRS-Komplex gehoren, 
dient Dazu zahlen die P- und T-Welle sowie hochfrequentere Storungen, 

25 die z. B. von der bioelektrischen Muskelaktivitat ausgehen. Der eingesetzte 
Filter BP ist linearphasig, nicht rekursiv und besitzt eine Bandpaficharakte- 
ristik mit den DurchlaBfrequenzen f g , = 18 Hz und fg2 ^27 Hz sowie den 
Sperrgrenzfrequenzen f sl = 2 Hz und fs2= 50 Hz. Die Filterordnung ist 
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FO = 26. Die Gruppenlaufzeit des BandpaBfilters BP entspricht dement- 
sprechend 13 Abtastwerten und muB bei der Bestimmung des Zeitpunktes 
des QRS-Komplexes beriicksichtigt werden. 

5 Die so erhaltenen Signalwerte x f (n) werden anschliefiend in einer Qua- 
drierstufe QS gem&B folgender Beziehung quadriert: 

I Xf(n)l 2 

10 Die so aus den urspriinglichen Signalwerten x(n) durch eine Art Betrags- 
bildung aufbereiteten Werte Xfq(n) werden nun einem Vergleicherkomplex 
HZ zugefuhrt, der diese Signalwerte mit einer adaptiv aus diesen ermittel- 
ten SchweUenwertfunktion K(n) vergleicht In Fig. 2 ist der die Ermittlung 
der Schwellenwertfunktion K(n) betreffende Verfahrenskomplex mit AS 
1 5 gekennzeichnet. Darin wird ein geeigneter Wert fur die Funktionskoeffizi- 
enten K(n) adaptiv aus den Signalwerten x fq (n) geschatzt. Dazu werden 
flieBend die BandpaB-gefilterten und quadrierten Signalwerte mit Hilfe ei- 
nes Gedachtnisfaktors 
A, k (0 < ^ k < 1) rekursiv gemittelt. 



K(n) = >, k K(n- 1) + (1- >. k )x fq (n)#c 
wobei c ein Konstante ist. 

25 Empirisch ergeben sich als geeignete Werte X k = 0,98 und c = 8. 



Die durch den GedachtmsfaktorA* vorgegebene Mittelungszeit bestimmt im 
wesentlichen die Adaptionsgeschwindigkeit dieser Schatzamg, wobei so- 
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wohl zu kiirze als auch zu lange Mittelungssegmente die Leistungsfahigkeit 
des Signalauswerteverfahrens beeintrachtigen konnen. 

Im Verfabxenskomplex HZ werden nun die Signalwerte x fq (n) mit der 
5 SchweUenwertfunktion K(n) - wie angesprochen - verglichen. Dabei wird 
festgestellt, in welcher Richtung die Signalwerte x fq (n) von der Schwellen- 
wertfunktion K(n) abweichen. Daraus wird eine Haufigkeitsanzahl D(n) 
innerhalb dieses definierten Segmentes N festgestellt, die die Anzahl oder 
Haufigkeit von Ereignissen reprasentiert, bei denen die Betrage der Signal- 
10 werte x fq (n) unterhalb der Schwellenwertfunktion liegen. D(n) kann auch 
auf rechengiinstige Weise rekursiv uber 



ermittelt werden. Ein kleiner Wert fur die Haufigkeitszahl D(n) spricht da- 
fur, daB der Betrag des EKG-Signals 4 dauerhaft oberhalb der Schwellen- 
wertfunktion K(n) liegt, was einen verlaBlichen Parameter fur das Vorlie- 
20 gen des QRS-Komplexes darstellt. 

Im Zuge des erfindungsgemaBen Verfahrens mufi also noch ein Grenzwert 
0 ermittelt werden, dessen Unterschreiten significant fur das Vorliegen 
eines QRS-Komplexes in dem definierten Segment des EKG-Signals 4 
25 spricht. 0(n) wird rekursiv aus D(n) berechnet mit 

0(n) = Xq 0(n - 1) + (1 - Xe) D(n) 



15 



D(n) = X D D(n - 1) + (1 - X D ) d(n) mit d(n) =. 




wobei ein Gedachtnisfaktor 0 < Xq < 1 verwendet wird. Fiir das am Ende 



dieser Beschreibung kurz wiedergegebene Zahlenbeispiel wurde dieser Ge- 
dachtnisfaktor zu Xe = 0,99 gewahlt. Liegt D(n) unterhalb des Grenzwertes 
0, ist ein QRS-Komplex erfafit worden, andernfalls nicht. 

Die Abprufiing, ob die vorstehende Bedingung erfiillt worden ist, findet in 
der Entscheiderstufe ES gemaB Fig. 2 statt. 
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Patentanspriiche 

1 . Signalauswerteverfahren zur Detektion von QRS-Komplexen in Elek- 
trokardiogramm-(EKG-)Signalen, gekennzeichnet durch folgende Ver- 
fahrensschritte: 

- Abtasten des EKG-Signals (4) und Umwandlung in diskrete, zeitlich 
aufeinanderfolgende Signalwerte (x(n)), 

- Vergleichen der Signalwerte (xf(n), x fq (n)) mit einer adaptiv aus die- 
sen ermittelten SchweUenwertfunktion (K(n)), 

- Feststellen einer Haufigkeitszahl (D(n)) innerhalb eines definierten 
Segmentes der aufeinanderfolgenden Signalwerte (xf(n), x fq (n)), mit 
der die Signalwerte (xf(n), x fq (n)) die Schwellenwertfunktion (K(n)) 
unterschreiten, und 

- Vergleichen der ermittelten Haufigkeitszahl (D(n)) mit einem defi- 
nierten Grenzwert (0), wobei ein Unterschreiten des Grenzwertes 
(0) signifikant fur das Vorliegen eines QRS-Komplexes (5, 6, 7) in 
dem definierten Segment des EKG-Signals (4) ist. 

2. Signalauswerteverfahren nach Anspruch 1, dadurch gekennzeichnet, 
daB die EKG-Signalwerte (x(n)) nach dem Abtasten einer HochpaBfilte- 
rung, vorzugsweise einer BandpaBfilterung (BP) unterworfen werden. 

3. Signalauswerteverfahren nach Anspruch 2, dadurch gekennzeichnet, 
daB die untere und obere DurchlaBgrenzfrequenz (f gU fg2) des BandpaB- 
filters (BP) bei ca. 18 Hz und ca. 27 Hz liegen. 

4. Signalauswerteverfahren nach einem der Anspruche 1 bis 3, dadurch 
gekennzeichnet, daB die Signalwerte (x(n)) einer Betragsbildung unter- 
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worfen werden. 

5. Signalauswerteverfahren nach Anspruch 3, dadurch gekennzeichnet, 
daC die Signalwerte (xf(n)) vor dem Vergleichen mit der Schwellenwert- 
fiinktion (K(n)) und Ermitteln der Haufigkeitszahl (D(n)) zur Betragsbil- 
dung quadriert werden. 

6. Signalauswerteverfahren nach einem der Ansprilche 1 bis 4, dadurch 
gekennzeichnet, daB der Wert der SchweUenwertfunktion (K(n)) adap- 
tiv aus einer fliefienden Mittelung der Signalwerte (x(n)), vorzugsweise 
der Bandpafi-gefilterten und quadrierten Signalwerte (x fq (n)), ttber eine 
durch den Gedachtnisfaktor (^ k ) vorgegebene Mittelungsperiode (P) be- 
stimmt wird. 

7. Signalauswerteverfahren nach einem der Anspruche 1 bis 5, dadurch 
gekennzeichnet, daB der fur einen QRS- Komplex (5, 6, 7) signifikante 
Grenzwert (0) der Haufigkeitszahl (D(n)) als adaptive Schwelle aus der 
Haufigkeitszahl (D(n)) selbst variabel eingestellt wird. 
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Zusammenfassung 

Ein Signalauswerteverfahren zur Detektion von QRS-Komplexen in Elek- 
trokardiogramm-Signalen weist folgende Verfahrensschritte auf: 

- Abtasten des EKG-Signals (4) und Umwandlung in diskrete, zeitlich 
aufeinanderfolgende Signalwerte (x(n), x fq (n)), 

- Vergleichen der Signalwerte (x(n), x fq (n)) mit einer adaptiv aus die- 
sen ermittelten SchweUenwertfiinktion (K(n)), 

- Feststellen einer Haufigkeitszahl (D(n)) innerhalb eines definierten 
Segmentes der aufeinanderfolgenden Signalwerte (x(n), x fq (n)), mit 
der die Betrage der Signalwerte (x(n), Xfq(n)) die Schwellenwert- 
funktion (K(n)) unterschreiten, und 

- Vergleichen der ermittelten Haufigkeitszahl (D(n)) mit einem defi- 
nierten Grenzwert (GW), wobei ein Unterschreiten des Grenzwerts 
(GW) signifikant fur das Vorliegen eines QRS-Komplexes (5, 6, 7) 
in dem definierten Segment des EKG-Signals (4) ist. 



-Fig. 2- 



